Comparison of the Legacy and Gold Snla Dataset Constraints on Dark Energy Models 
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We have performed a comparative analysis of three recent and reliable Snla datasets available 
in the literature: the Full Gold (FG) dataset (157 data points < z < 1.75), a Truncated Gold 
(TG) dataset (140 data points < z < 1) and the most recent Supernova Legacy Survey (SNLS) 
dataset (115 data points < z < 1). We have analyzed and compared the likelihood of cosmological 
constant and dynamical dark energy parametrizations allowing for crossing of the phantom divide 
line (PDL). We find that even though the constraints obtained using the three datasets are consistent 
with each other at the 95% confidence level, the latest (SNLS) dataset shows distinct trends which 
are not shared by the Gold datasets. We find that the best fit dynamical w(z) obtained from the 
SNLS dataset does not cross the PDL w = — 1 and remains above and close to the w = — 1 line 
for the whole redshift range < z < 1 showing no evidence for phantom behavior. The LCDM 
parameter values (wo = — 1, Wi = 0) almost coincide with the best fit parameters of the dynamical 
w(z) parametrizations. In contrast, the best fit dynamical w(z) obtained from the Gold datasets 
(FG and TG) clearly crosses the PDL and departs significantly from the PDL w = — 1 line while 
the LCDM parameter values are about 2a away from the best fit w(z). In addition, the (fiom,£iA) 
parameters in a LCDM parametrization without a flat prior, fit by the SNLS dataset, favor the 
minimal fiat LCDM concordance model. The corresponding fit with the Gold datasets mildly favors 
a closed universe and the flat LCDM parameter values are la — 2a away from the best fit (fiom, ^a). 

PACS numbers: 98.80.Es,98.65.Dx,98.62.Sb 



I. INTRODUCTION 

Current cosmological observations show strong evi- 
dence that we live in a spatially flat universe [l[ with 
low matter density Q that is currently undergoing ac- 
celerated cosmic expansion 0, 0, @- The most direct 
indication for the current accelerating expansion comes 
from the accumulating type la supernovae (Snla) data 
0, Q which provide a detailed form of the recent expan- 
sion history of the universe. 

This accelerating expansion has been attributed to a 
dark energy 6] component with negative pressure which 
can induce repulsive gravity and thus cause accelerated 
expansion. 

The simplest and most obvious candidate for this dark 
energy is the cosmological constant A _7j with equation 
of state w = p/ p = — 1. This model however raises theo- 
retical problems related to the fine tuned value required 
for the cosmological constant Q- These difficulties have 
lead to a large variety of proposed models where the dark 
energy component evolves with time |8j usually due to an 
evolving scalar field (quintessence) which may be mini- 
mally ||| or non-minimally [t| coupled to gravity. The 
main prediction of the dynamical models is the evolution 
of the dark energy density parameter Clx(z). Combining 
this prediction with the prior assumption for the matter 
density parameter f2om, the predicted expansion history 
H(z) is obtained as 



as 
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If the dark energy can be described as an ideal 
fluid with conserved energy momentum tensor T^ u = 
diag(p,p,p,p) then the above parameter w(z) is identi- 
cal with the equation of state parameter of dark energy 



w(z) 
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H(z) 2 = H^[n 0m {i + z) 3 + n x (z)} 



(i.i) 



The dark energy density parameter is usually expressed 



Independently of its physical origin, the parameter w(z) 
is an observable derived from H(z) (with prior knowledge 
of f^om) and is usually used to compare theoretical model 
predictions with observations. 

Most evolution behaviors of w(z) can be reproduced by 
assuming appropriate scalar field quintessence potentials. 
If however w(z) were observationally found to cross the 
phantom divide line (PDL) w = — 1 then all minimally 
coupled single scalar field models would be ruled out as 
dark energy candidates IT^. (this includes phantom [l4| 
and k — essence models |15|). This would leave only mod- 
els based on extended gravity theories 0, 0] and com- 
binations of multiple fields (la, Il9j (quintessence + phan- 
tom) as dark energy candidates. It is therefore important 
to utilize the available Snla datasets to place constraints 
on w{z) and determine the likelihood of having a w{z) 
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that crosses the PDL, as done in Ref.po| using the Gold 
dataset following a robust procedure. This paper con- 
firms the evidence for crossing the PDL with w < — 1 at 
present and also shows some hint for oscillating w{z) at 
best fit in agreement with Ref.|ll|. which were the the 
first to point this trend. 

The two most reliable and robust Snla datasets ex- 
isting at present are the Gold dataset [4( and the Su- 
pernova Legacy Survey (SNLS) [f| dataset. The Gold 
dataset compiled by Riess et. al. is a set of supernova 
data from various sources analyzed in a consistent and 
robust manner with reduced calibration errors arising 
from systematics. It contains 143 points from previously 
published data plus 14 points with z > 1 discovered re- 
cently with the HST. The SNLS is a 5-year survey of 
Snla with z < 1. It has recently Q released the first 
year dataset. The SNLS has adopted a more efficient 
Snla search strategy involving a 'rolling search' mode 
where a given field is observed every third or fourth night 
using a single imaging instrument thus reducing pho- 
tometric systematic uncertainties. The published first 
year SNLS dataset consists of 44 previously published 
nearby Snla with 0.015 < z < 0.125 plus 73 distant Snla 
(0.15 < z < 1) discovered by SNLS two of which are out- 
liers and are not used in the analysis. The fact that in 
the two datasets a set of low-z Snla is common to both 
samples could only lead to minor common systematics 
due to low redshift. 

The above observations provide the apparent mag- 
nitude m(z) of the supernovae at peak brightness af- 
ter implementing correction for galactic extinction, K- 
correction and light curve width-luminosity correction. 
The resulting apparent magnitude m(z) is related to the 
luminosity distance cIl(z) through 

m th (z) = M(M,H ) + 5log 10 (D L (z)) (1.5) 

where in a flat cosmological model 

Dl{z) = (1 + z) f dz'-—^ T (1.6) 

J H(z';ai, ...,a n ) 

is the Hubble free luminosity distance (HqcLl/c), 
Oi, a n are theoretical model parameters and M is the 
magnitude zero point offset and depends on the absolute 
magnitude M and on the present Hubble parameter Hq 
as 

c Ff- 1 

M = M + 5logi Q (—Q-) + 25 = 
Mpc 

= M — 5log w h + 42.38 (1.7) 

The parameter M is the absolute magnitude which is 
assumed to be constant after the above mentioned cor- 
rections have been implemented in m(z). 

The data points of the Gold dataset are given after 
the corrections have been implemented, in terms of the 
distance modulus 

^ bs (zi) = m^ s (zi)-M (1.8) 

The SNLS dataset however also presents for each point, 
the stretch factor s used to calibrate the absolute magni- 
tude and the rest frame color parameter c which mainly 



measures host galaxy extinction by dust. Thus, the dis- 
tance modulus in this case depends apart from the abso- 
lute magnitude M , on two additional parameters a and 
(3 defined from 

»Ts LS = m S Z LS { Zl ) - M + a( Sl - 1) - (5 Cl (1.9) 

which arc fit along with the theoretical parameters using 
a recursive procedure discussed below. 

The theoretical model parameters are determined by 
minimizing the quantity 

X (Oi.-.ftn) = „2 ,„2 ,„2 ( L1 °) 

i=l f i int v i 

where a 2 i; a 2 nt and a 2 , are the errors due to flux uncer- 
tainties, intrinsic dispersion of Snla absolute magnitude 
and peculiar velocity dispersion respectively. These er- 
rors are assumed to be gaussian and uncorrelated. The 
theoretical distance modulus is defined as 

Hth(zi) = m th (zi) - M = Uog w (D L {z)) + ^ (1-H) 
where 

Mo = 42.38 - 5log 10 h (1.12) 

and n bs is given by ljl.8|l and H1.9|) for the Gold and 
SNLS datasets respectively. 

The steps we followed for the minimization of (|1.10|) 
for the Gold and SNLS datasets are described in detail 
in the Appendix. The validity of our analysis has been 
verified by comparing the part of our results that overlaps 
with the results of the original Refs [4J, |5j of the Legacy 
and Gold datasets. 



II. COMPARATIVE ANALYSIS 

We will consider four representative H(z) parametriza- 
tions and minimize the % 2 of eq. Hl.lOfl with respect to 
model parameters. We compare the best fit parametriza- 
tions obtained with three datasets 

• The full SNLS dataset with 115 datapoints (exclud- 
ing two outliers) and z < 1. 

• The Full Gold dataset (FG) with 157 datapoints 
and < z < 1.7. 

• A Truncated version of the Gold dataset (TG) with 
140 datapoints and z < 1 which can be compared 
in a more direct way with SNLS. 

The four fitted parametrizations include the general 
LCDM without a flat prior: 

H{z) 2 = H 2 [Q 0m {l + zf + n A + 

+ (i - n 0m - O0(i + z) 2 } (2.i) 

w{z) = -1 (2.2) 

and three dynamical dark energy parametrizations with 
two free parameters which allow for crossing of the PDL: 
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fi M = 0.76 ± 0.24 
n A = 1.40 + 0.35 
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FIG. 1: The 68% and 95% confidence region ellipses in the f2nm — plane based on parametrization 12. II . The three plots 
correspond to the three datasets discussed in the text (SNLS, TG and FG). Notice that flat LCDM is more favored by the 
SNLS dataset than by the Gold datasets. 



• Parametrization A: 

w(z) = wq + w\ z (2-3) 
H 2 (z) = H 2 [Q Qrn (l + zf + (2.4) 
+ (l-^ m )(l + z) 3(1+lU0 -™ l) e 3 " ,lZ ] 



• Parametrization B: 



w{z) = Wq + W\- 



l + z 

H 2 {z) = Hl[VL 0m {l + zf + 
+ (1 - fi m)(l + z) 3( - 1+W0+w ^e 3wi[1/il+z) - 1] ] (2.6) 

• Parametrization C: 

, ai + 3(Q m - 1) - 2aiz - a 2 (-2 + 2z + z 2 ) 

w ( z > = 571 o Z Zo Z T\ ( 2J ) 

3(1 - S2oro + a i z + 2a 2 z + a 2 z z ) 

H 2 (z) = H 2 [n 0m (l + zf + ai(l + z) + 02(1 + z) 2 + 

+ (1 - O 0m - 01 - a 2 )} (2.8) 

The motivation behind parametrization C is to mimic 
a two-component DE model. Alternatively, it could be 
viewed as a power law expansion in the scale factor de- 
pendence of the DE energy density. In analyzing the 
general LCDM of ea. (|2.1() we used 



(2.5) 



Dl(z) 



(1 + *) 



yOom + ^ox — 1 



Sin[y/tt 0m + ^ox - 1 / dz 



(2.9) 



H(z)' 



instead of eci. (|1.6|) which is only suitable for flat models. 

In Fig. 1 we show the confidence region ellipses in the 
Oom — plane based on parametrization (|2.1() . The 
three plots correspond to the three datasets discussed 
above (SNLS, TG and FG). The following comments can 
be made on these plots: 



• The major axes of the elliptic contours are approx- 
imately parallel in the three plots. This effect [2 lj is 
due to the degeneracy of the fitted Dl(z; flo m , ^a) 
with respect to certain linear combinations of the 
parameters Qo m — JIa in the redshift range of inter- 
est. For example, choosing a representative red- 
shift z — 0.6 it is easy to show that the value 
Dl(z = 0.6) = 0.84 is obtained by all combinations 
of flom — that satisfy Qo m — 0.80fiA = —0.38. 
The direction of this (approximate) degeneracy line 
is determined by the H(z) parametrization and the 
redshift range considered but the actual location of 
the line is determined by the data. 

• The two versions of the Gold dataset favor a closed 
universe instead of a flat universe = 2.16 ± 
0.59, 0,[ o f = 1.44±0.44). This trend is not realized 
by the SNLS dataset which gives nf ^ LS = 1.07 ± 
0.52. 

• The point corresponding to SCDM (f2om,^A) = 
(1,0) is ruled out by all datasets at a confidence 
level more than 10cr. 

• The values of /i™ m that minimize the x 2 (^0m, ^a) 
of eq. i|1.10fl (with H(z) given by obtained 
by all three datasets are consistent with each other. 
We find ^ NLS = 43.15 ± 0.05, /i^ G = 43.30 ± 
0.05 and ^ G = 43.32 ± 0.05. This alleviates the 
possible discrepancy between high and low redshift 
data discussed recently in Ref. (21| . 

• If we use a prior constraint of flatness ilom + ^A = 1 
thus restricting on the corresponding dotted line of 
Fig. 1 and using the parametrization 



H(zf = H 2 [n 0m (i + zf + (1 - n 0m )} 



(2.10) 
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FIG. 2: The 68% and 95% \ 2 confidence contours of parametrizations A and B assuming a flat prior using the datasets SNLS, 
TG and FG. A prior of Qom ~ 0.24 has been used. The dashed lines intersect at the parameter values of flat LCDM. Notice 
that for the SNLS, flat LCDM almost coincides with the best fit for both parametrizations. 



we find minimizing x 2 (^0m) of eq (|1 . 1 C)|) 

n^ LS = 0.26 ±0.04 (2.11) 
fi£G = o.30 ±0.05 (2.12) 
n££ = 0.31 ±0.04 (2.13) 

The values of ^Iq^ ls and fig® are practically iden- 
tical with the corresponding in the original Refs 
0,O| where the data were first published. 

This, along other similar tests, confirms the validity 
of our analysis. 

Even though the cosmological constant with equation 
of state parameter w = — 1 is the simplest form of dark 
energy consistent with the data, the possibility of evolv- 
ing dark energy models with non-constant w(z) remains 
a viable alternative which may even provide better fits to 
the data than LCDM. To address this issue we considered 
the three parametrizations (A, B, C) of eqs l|2.3|l - (|2.8|) 



and assuming flatness we constrained their parameters 
using the three datasets. 

Parametrizations A and B reduce to w(z) = wq in the 
special case when w does not evolve with time (wi = 0). 

With this prior we construct Table I showing the best 
fit wq value and la errors for each dataset and flo m — 
0.24 — 0.3. Clearly all datasets are consistent with each 
other ruling out models with w > 1/3 at more than a 
10<7 confidence level. 

Extending the analysis to the full parameter space we 
construct x 2 confidence contours in Figs. 2 and 3 assum- 
ing a flat prior. 

In order to investigate the dependence of our results on 
the prior of f^om we do not marginalize over it. Instead we 
present the wo — w\ confidence contours of parametriza- 
tions A and B for the priors Qo m = 0.24 (Fig. 2) and 
Uom = 0.3 (Fig. 3). The point corresponding to LCDM 
is also shown in these figures. For parametrization C we 
have not presented the confidence contours because H{z) 
becomes complex for regions of parameter space overlap- 
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FIG. 3: Same as Fig. 2 but with a prior Qo m = 0.3. 



ing with the 95% confidence contours. 

The following comments can be made regarding Figs. 
2 and 3. 

• The minimal LCDM model (wq = —1, w\ = 0) 
appears to be close to the 95% confidence contour 
in the analyses based on the TG and FG datasets. 
For the analysis based on the SNLS however, the 
flat LCDM model is well within the 68% contour 
and in fact for f2 Um = 0.24 it is almost identical 



with the best fit paramctrization in both the A and 
B parametrization cases! Thus LCDM appears to 
have significantly gained in likelihood compared to 
dynamical dark energy models in the context of the 
new SNLS dataset. 



TABLE I: The best fit wo parameter values for each dataset 
assuming priors of flatness and wi — 0. 



Matter Density Qom 


SNLS 




TG 


FG 


0.24 w = 


= -0.95 ± 


0.09 


too = -0.89 ± O.fO 


w = -0.86 ± 0.09 


0.30 w = 


= -l.fl ± 


O.ff 


w = -1.04 ±0.12 


w = -1.02 ±0.11 



6 



• The best fit w(z) parametrizations in the context 
of the FG and the TG datasets, not only are far 
from the LCDM point (wq = — 1, wi = 0) (see 
Figs. 2 and 3) but they also clearly cross the PDL 
line. This is demonstrated in Fig. 4 where the 
best fit w{z) are plotted for each dataset and each 
parametrization with a prior of Oo m = 0.24. 

This crossing of the PDL is not realized for the 
best fit A, B and C parametrizations in the context 
of the SNLS dataset. Several authors^! E2 
have been recently motivated by the high likeli- 
hood of the PDL crossing indicated by the Gold 
datasets |23j, |24| to explore theoretical models that 
predict such crossing. It has been shown that this 
task is not trivial and can not be achieved by a 
single minimally coupled field 0. 

Here we show however that such PDL crossing is 
not favored by the new SNLS dataset and therefore 
the motivation for the above papers is weakened. 
Indeed, phantom[14| dynamical dark energy models 
with wq < — 1 are not favored by the SNLS dataset 
in contrast with the Gold dataset that favored such 
models (see Figs 2 and 3 and Ref. [ME^L 

• Even though the best fit parameter values are rel- 
atively insensitive to the value of ilom m the range 
0.2 < Qom < 0.3 the errors are more sensitive in- 
creasing with Qom- Our results for the best fit pa- 
rameters of the parametrizations A, B and C are 
shown in the legends of Figs 2 and 3. For both 
JJom priors the errors of w\ are much larger com- 
pared to those of wq. As discussed above, this is a 
generic feature of the parametrization used and is 
not related to the particular datasets. 

III. DISCUSSION-CONCLUSION 

We have performed a comparative analysis of the three 
most recent and reliable Snla datasets available in the 
literature: the Full Gold (FG) dataset (157 data points 
< z < 1.7), the Truncated Gold (TG) dataset (140 data 
points < z < 1) and the most recent SNLS dataset (115 
data points < z < 1). Our analysis is an extension of 
our earlier analyses which had focused |25j on the FG 
and earlier datasets. We have used representative 
dark energy parametrizations to examine the consistency 
among the three datasets in constraining the correspond- 
ing parameter values. We have found that even though 
the constraints obtained using the three datasets are con- 
sistent with each other at the 95% confidence level, the 
latest (SNLS) dataset shows distinct trends which are not 
shared by the other (earlier) datasets. The most charac- 
teristic of these trends are the following: 

• The best fit dynamical w(z) obtained from the 
SNLS dataset does not cross the PDL w = — 1 and 
remains close to the w = — 1 line for the whole 



redshift range < z < 1. The LCDM parameters 
(wq = — 1, w\ = 0) almost coincide with the best 
fit parameters of the dynamical w(z) parametriza- 
tions. Thus, the data do not seem to require and 
utilize the additional dynamical parameters offered 
to them. This is an interesting new feature of the 
data which favors the minimal LCDM model. In 
contrast, the best fit dynamical w(z) obtained from 
the Gold datasets (FG and TG) clearly crosses 
the PDL and departs significantly from the PDL 
w = — 1 line (see Fig. 4). According to these 
datasets the minimal LCDM is consistent but is 
not favored. It is about 2a away from the best fit 
w{z) which crosses the PDL. 

• The best fit (SIojth^a) parameters in a LCDM 
parametrization without a flat prior show inter- 
esting differences between the Gold and the SNLS 
datasets. In particular, the SNLS favors a flat uni- 
verse much more than the Gold datasets. How- 
ever, all three datasets remain consistent with flat 
LCDM at the 95% confidence level while SCDM is 
excluded by all datasets at more than 10<r. 

The above mild trend differences between the Gold and 
the SNLS datasets can be summarized by stating that the 
SNLS hints towards the minimal flat LCDM concordance 
model more than the Gold datasets. It is an exciting 
prospect to see whether this trend will continue and get 
verified by upcoming future Snla observations. 



IV. APPENDIX 

Here we give some details of our analysis for both the 
Gold and the Legacy datasets. The full numerical analy- 
sis was performed using Mathematica and it is available 
at |htt p: / /leandros. physics. uoi.gr /snls.htm 

We first describe the method used for the Gold dataset 
analysis. From eqs. (|1.10|) . I|l.llfl and (|1.12|) we find 

2/ \ (Hobs i - 5log w D L (zi; au a n ) - ^ ) 2 

X [flu -,a n ) = 2 

i=i a * 

(4.1) 

where 

of =°li + aia + °l i ( 4 -2) 

is the total error published for the Gold dataset. 

The parameter /iq is a nuisance parameter but it is 
independent of the data points and the dataset. This 
expected independence can be used as consistency test 
of the data |2l|. The minimization with respect to /io 
can be made trivially by expanding the \ 2 of equation 
(|1.10[) with respect to fi as 

X 2 (au-,a n ) = A ~ 2fi B + v 2 C (4.3) 
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FIG. 4: The best fit w(z) are plotted for each dataset (SNLS, TG and FG) and each parametrization (A, B and C) with a 
prior of Slo m = 0.24. The thick solid line represents the best-fit and the light grey contour represents the la confidence level 
around the best-fit. The dased horizontal line represents LCDM. Notice that PDL crossing at best fit occurs only for the TG 
and the FG (Gold) datasets while the best fits of the SNLS dataset almost coincide with LCDM. 



where Equation (|4.3|l has a minimum for fiQ = B/C at 
A( ail ..,a rl ) - ^ ("w(^)-"W^o = 0,ai,.,a»)) 2 jffa, a n ) = A(oi, a n ) - a ^ (4.5) 

1=1 a m obs (zi) 

N Thus instead of minimizing x (MO) a i, •••> ^n) we can mini- 
D/ s (m obs (zi) - m t h{zi\no = 0,oi,.. 

i a n)) mize x 2 ( a i: a n) which is independent of /io- Obviously 

<J m obs (z i ) Xmin = Xmin an( i this ^ s * ne approach used. Alterna- 

N tively we could have performed a uniform marginalization 

Y,^— (4-4) 

i=l m obB (zi) 



c 



8 



over the nuisance parameter /j,q thus obtaining [TTL fl3L l2^ | 



-2/ \ a, \ B(a 1 , a n ) 2 

X [at, ..., a n ) = A{ai, ...,a n ) — hln(G/27r) 

(4.6) 

to be minimized with respect to ai,...,a n . In our Gold 
dataset analysis we consider the x 2 (ai, a n ) of equation 
(14. 5|) which is already minimized with respect to /io- If 
we marginalized over all values of Hq, as in Ref. £|, that 
would just add a constant (see ea. l^.tif) ') and would not 
change the results. The minimization of l|4.1|l was made 
using the FindMinimum command of Mathematica. 

Our analysis of the SNLS dataset proceeded in a some- 
what different manner following Ref. p| . Using equation 
(11.91) in 11.10|l we constructed x 2 as 



N 

E 



(Hob 



X 2 (a,f3, M + ^ Q ,ai,...,a n ) = 
5log 10 D L (zi; ai, a n ) - /i ) 2 



(4.7) 



For the distance modulus error we have 

^^^L + A^ + Ao (4.8) 

Each one of the <j m i, a s i and c c i has been published 
in Ref. ||. The velocity dispersion error er 2 assuming a 
peculiar velocity dispersion of 300km/ sec may be written 
as 



5 10- 



1 



1 



(4.9) 



value <7i n t — 0.15 and then updated with the following 
three step procedure 

1 . Fix the values of a and (3 in cr 2 i and minimize the 
X 2 of eq. (|4.8|) with a int — 0.15. If this fixing is 
not made, a bias is introduced towards increasing 
errors during minimization. 

2. Change the value of Ui n t to obtain x 2 = 1- 

3. Use the new value of <Ji n t and minimize again keep- 
ing the values of a and (3 in er 2 i fixed. 

This procedure, with no marginalization over M + [j,q, a, 
(3 as described in page 10 of Ref. [5j, leads to the best fit 
values of M + fio, a , (3, ai,...,a n . 

The errors are evaluated using the covariance matrix of 
the fitted parameters [27J and the errors on the equation 
of state w(z;pi) are given by 



The intrinsic dispersion error <jf nt is initially set to a 



lLr\ ft fli i"\ i~\ 

^=E(^) C '«+ 2 E E ^9f? Cij (4 ' 10) 

i— 1 1=1 J 

where pi are the cosmological parameters and Cij the 
covariance matrix [2^] . 
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